Ring polymers in melts and solutions: scaling and crossover 
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We propose a simple mean-field theory for the structure of ring polymer melts. By combining the 
notion of topological volume fraction and a classical van der Waals theory of fluids, we take into 
account many body effects of topological origin in dense systems. We predict that although the 
compact statistics with the Flory exponent v = 1/3 is realized for very long chains, most practical 
cases fall into the crossover regime with the apparent exponent v = 2/5 during which the system 
evolves toward a topological dense-packed limit. 
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Statistics of melts and concentrated solutions of 
ring polymers is a longstanding problem in polymer 
physics p]4l4|. Unlike linear chain systems, their physi- 
cal properties crucially depends on the preparation his- 
tory during which the topology of the system is frozen. 
Then, the non-crossing requirement creates topological 
constraints which impose nontrivial restrictions on the 
phase space of the system, hence, have a strong influence 
even in statics. 

In the linear polymer counterpart, a well-known Flory 
theorem states that the chain conformation is Gaussian 
characterized by the Flory exponent v = 1/2 due to the 
screening of excluded- volume interactions |15Ml7l |. This 
simple (but surprising) result, however, no longer holds 
for the ring polymer melt. The most basic question arises 
in the melt of ring polymers free from any mutual-linking 
and self-knotting. There have been several experimen- 
tal @, 0| and numerical studies [8l-[l4j in this direction, 
but the clear-cut conclusion has not been attained yet. 

In their seminal paper, Cates and Deutsch (C-D) ar- 
gued that the unconcatenated rings in the melt may 
have statistics intermediate between those of collapsed 
(y = 1/3) and Gaussian (v = 1/2) chains [l[. Specif- 
ically, they proposed a conjecture on the scaling expo- 
nent v = 2/5 based on Flory-type mean-field theory. 
While this leading theoretical guide seems to be sup- 
ported by following numerical simulations some 
of more recent observations claim the collapsed statistics 
v = 1/3 for sufficiently long rings [TH, [lj] . The latter 
result was hypothesized by Khokhlov and Nechaev based 
on the analogy with lattice animals [|| . Closely related to 
this is the crumpled globule (CG) model [1,01 which was 
originally proposed by Grosbcrg, ct. al. as a long-lived 
kinetic intermediate on the pathway of the linear chain 
collapse [H , and then, hypothesized as a large scale DNA 
organization in interphase chromosomes. This suggests 
an interesting link between polymer topology and biolog- 
ically important problems such as the formation of the 
chromosome territories UM ■ 
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It is important to keep in mind that the topological ef- 
fect manifests itself in the scale larger than some charac- 
teristic length £i. Individual rings in concentrated solu- 
tions of small molecular weight N < gi = (Ci/a) 2 ^ 1 / 4 (N 
is the number of monomers in each ring, a is the monomer 
size and 4> is the monomeric volume fraction) thus show 
Gaussian behaviors with the size Rq — ^(N/g) 1 / 2 ~ 
aN 1 / 2 ^ 1 / 8 where £ ~ ag 3 ^ 5 ~ a</>~ 3 / 4 is the correla- 
tion length of concentration fluctuation. In the present 
paper, we introduce the notion of the topological volume 
fraction and construct a mean-filed theory for the con- 
centrated solution of noncatenated long ring polymers. 
Requiring the theory to be compatible with the above 
feature associated with the length scale £i , we show that 
the scaling exponent for the long chain limit is given by 
v = 1/3. However, extremely long chains N > N* = Cg\ 
with a large numerical factor C ~ 10 are required for this 
asymptotic to be reached, thus, most practical cases fall 
into a broad crossover region (gi < N < N*) where the 
apparent exponent is given by v — 2/5. At various stage 
of the paper, we will look into the physics behind the C-D 
theory in the light of the present argument. We also sug- 
gest a connection to the CG model through a topological 
blob which adjusts its own size in an intriguing way in 
the crossover region. This view provides a natural bridge 
between otherwise conflicting two lines of previous con- 
jectures (C-D and CG). 

Mean field theory — A rigorous analysis of the topo- 
logical effects may rely on the mathematical knot invari- 
ance, which, however, seems to be yet formidable. We 
instead seek for a physically motivated coarse-grained 
phenomenological description. A successful prototype of 
such an approach is already found in the dynamics and 
rhcology of high molecular weight linear polymer melts in 
which the tube model provides a bridge between molecu- 
lar level topological constraints and macroscopic material 
properties (l5l . fl6| . In the present problem, we also pos- 
tulate an intrinsic length scale £i which is an analogue 
of the tube diameter in the linear chain solutions, below 
which the topological effects are irrelevant. At this stage, 
what is needed is an element to treat large scale static 
behaviors of unconcatenated rings, just as the reptation 
theory provides a basic framework to treat the dynamical 
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properties of entangled linear chain solutions beyond the 
tube size. 

A key observation from previous studies is that the 
topological constraints could be effectively represented 
as excluded-volume effects. Notable topics in this line in- 
clude an entropic repulsion between untangled ri ngs I19B . 
a topological swelling of isolated random knots [20, [21j 
and the anomalous bond correlation function of planar 
rings [22j]. It is expected that the above conjecture, first 
proposed by des Cloizeaux [23|, would be equally useful 
for the present problem, too. 

Let us consider a concentrated solution of ring poly- 
mers. The equilibrium spatial size R (such as a radius of 
gyration) of each ring is a function of N and the monomer 
volume fraction </>. The number of total monomers in the 
region of volume ~ R 3 of single ring is ~ R 3 (f>/a 3 , thus, 
the number of rings there is 

<» 

= a 3 N/R 3 , it can be written as 



Using a "self density" (j) 

N R = <j>/<f> s . 

The topological contribution to the free energy consists 
of terms arising from (i) noncatenation constrains among 
different rings and (ii) intra-ring effects associated with 
self- knotting 

R R-inter "h Rintra (2) 

In C-D theory, the first term is evaluated by assigning 
~ k B T to each of Nr neighboring rings, leading to the es- 
timate Pinter /ksT — {R? ' <j>) / '{a? ' N) . One may notice that 
this amounts to the evaluation of binary interactions. In 
concentrated solutions of long rings, however, the many- 
body correlation effects should become progressively im- 
portant, which may be taken into account in line with 
the excluded-volume analogy as followings. The effective 
excluded-volume vr of the ring with the spatial size R 
scales its volume 

v R ~ R 3 Y (3) 

where a dimensionless factor Y (independent of N) shall 
be fixed later to be consistent with the presence of the 
length scale £1 mentioned earlier. To evaluate the re- 
sultant repulsive interaction, we adopt a classical van 
der Waals theory of fluids whose free energy density 
for a one component fluid with the volume fraction <p 
and the excluded volume v is given by f '((j))/ '(k B T) = 
[4>/v\ In [0/(1 — 4>)] — e<f) 2 . In our case, this may be trans- 
formed as j ' inter /(k B T) ~ -[<j> R /v R ]ln[(l - 4>r)} where 
we set e = (at hernial), and as usual for the polymer 
solution theory, the ideal gas term is irrelevant here [l7| , 
and 



vrN r 
R 3 



= N R Y 



(4) 



is the ring's "volume fraction" of topological origin. The 
free energy per ring P in ter = f inter x R 3 /N R is thus 



k B T 



ln(l 



(5) 



An inspection of eq. ([5]) indicates that the larger size R 
costs free energy due to the noncatenation constraint. 
Thus, this topology effect leads to the squeezing of the 
ring toward a globular state, which should be negotiated 
with the unknotting constraint. A scaling analysis sug- 
gests that this free energy cost for the confinement should 
be written as 
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where R = ^(N/g) 1 / 2 = aN 1 / 2 ^ 1 ^ is the size of un- 
perturbed Gaussian ring (more precisely, the Gaussian 
ring of the blobs of size £ = ag 3 ^ 5 = a4>~ 3 ^ 4 ), and use 
has been made of eqs. (U} and (fj| in the last equality. 
One then requires Rintra to be extensive (~ TV) with 4>r 
given, which leads to 5 — 6 [25] ; 



Rtntra{4>R] N ) 
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(7) 



Note that this differs from a naive guess 8 = 2 adopted in 
the C-D theory which, however, is no longer correct for 
non-ideal chains confined in closed cavity (see ref. [24| 
for a general discussion on it). Indeed, eq. ([7]) is nat- 
urally expected in the context of the excluded-volume 
analogy [PA |2~H and applies when rings feel squeezing 
R < Rq O N > g\. It is noted that while Pinter solely 
depends on the topological volume fraction <f>R, Rintra 
linearly increases with the chain length N at a given <f>R. 
This reflects different aspects of the constraints between 
nonconcatenation and unknotting. 




FIG. 1: Plots of free energy as a function of <f>u for short ring 
BN/g-i = 0.4 (left) and long ring BN/gi = 10 (right). 

a. Short scale behavior Figure [1] shows typical free 
energy profiles for short and long rings. For small N, the 
free energy minimum is attained at low <j)R < 1. One 
may find a situation Pi n ter{(t>B,) - F intra {4> R ; gi) ~ k B T 
at 4>r. This indicates that, in short length scale r < 



= ar,}' 2 ^/* 



ii = a9i 



the topological effect is not effective, 



and the ideal ring behavior is observed. From the above 
condition, we find 



Y^ 9 - 1/2 r b/s S 1/2 (0fl) 



(8) 
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where B(4>]i) ~ <fi R ~ 0.5. We shall shortly argue that 
this numerical factor B is crucial for the characterization 
of the crossover regime. The confinement free energy 
eq. ([7]) can now be rewritten as 

FjntrajfoN) 2 TV 2 TV 

gi 92 

where g 2 = gi/B. 

b. Free energy minimization Now let us seek for the 
spatial structure of rings in large scales TV > g\. Here, 
the energy scale is much larger than the thermal energy 
(see Fig. []]), thus, the equilibrium size would be readily 
determined through the free energy minimization with 
respect to <pR. The result of numerical solution shown 
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FIG. 2: Normalized plot of the ring size R/(,2 as a function 
of the chain length N/g2 (double logarithmic scale). 

in Fig. [5] is summarized as follows; (i) sufficiently long 
rings (TV > TV*) obey a compact statistics with v = 1/3, 
(ii) rings with intermediate length (gi < TV < TV*) are 
well characterized by the effective exponent 2/5. There- 
fore, one may say that the statistics of concentrated ring 
polymers has well-defined two regimes (ideal and com- 
pact statistics for short TV < gi and long TV > TV* rings, 
respectively) which are separated by a wide crossover re- 
gion spanning over about one order of magnitude. To 
get an insight into this crossover, we note that from 
cqs. dTJ, CP, ©,©,©, a mean-field solution for TV > g 1 
can be generally written in the form 

where £ 2 = acj)" 1 ^ 'g^ 2 = £iB _1 / 2 and cj>^ denotes the 
equilibrium value which depends solely on N/g 2 . This 
indicates that upon renormalization a large scale spa- 
tial organization of individual rings is formally charac- 
terized by a space-filling curve with fractal dimension 




FIG. 3: Dependence of the self density <j> a and the number of 
neighboring rings TVjj on the chain length N (semi-logarithmic 

plot). 4> a decreases from <fi a ~ (j^^g^ 1 ^ 2 at TV = g\ toward 
4>a R ~ 1 ' — </> 3//8 <7^" 1//2 at TV > TV*, whereas Nr increases from 
TV H ~ <j) b/s g\ /2 at TV = gi toward TVjf n=1) ~ (p 5/s gl /2 at 
TV > TV*. 



v 1 = 3. Unlike usual cases, however, a renormaliza- 
tion factor, which is invoked in an earlier study as topo- 
logical blob £\)(N) = agl (N)4>~^/ & is not constant, 
but rather depends on TV, i.e., starting from g§ = g\ 
at TV = gi, it slowly (logarithmically) increases toward 
g\\ = .92 at TV = TV*, then saturates. It is this slow evo- 
lution toward a topological dense-packed limit (cf>R — > 1) 
that renders a small excess correction to the exponent in 
the large crossover region. Such a feature is demonstrated 
in Fig. [3] where the self-density <f) s and Nr (roughly corre- 
sponding to the number of neighboring rings) are plotted 
against the chain length. A point TV = TV* at which these 
observables saturate corresponds to the onset of the com- 
pact statistics, i.e., c/)r — > 1 at TV > g 2 which leads to the 
onset length TV* = Gg x with C ~ 10 (Fig. ©. To get a 
number, let us adopt an entanglement length as a mea- 
sure of <7i, and substitute a typical value ~ 100. This 
yields TV* - 1000 

Several remarks are now in order. As a mean-field 
theory, if appropriately constructed, is expected to pro- 
vide an accurate description on concentrated solutions, 
it would be of value to get another perspective on the 
structure of free energies. First, eq. © for F inter can be 
expanded in a virial series Fi nter /ksT — 0fi + </>lj/2 + - • • . 
Retaining only a lowest order term (second virial approx- 
imation), one finds Fi nter /kBT ~ R 3 (f>Y/(a 3 N) which 
corresponds to the free energy adopted in the C-D the- 
ory aside from a factor Y [27|. Our analysis suggests, 
however, that higher order terms are essential. Sec- 
ond, a view of the compact ring of size R (< Rq) as 
a dense piling of topological blob implies ^ ~ S.ig^/g) 1 ^ 2 
and a 3 g^/£ 3 = <j> s . One thus finds ^ ~ a0 1 / 4 (/)J 1 ~ 
a(t)~ 3/4 N R and ~ _5 / 4 TV|.. Assigning ~ k B T per 
topological strand, we obtain Fi ntra / {ksT) ~ N/g^ ~ 
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R 3 /£ 3 , which coincides with the scaling derivation of 
Fintra (cq. ©) [24| . Moreover, the deduced relations 
^ ~ £Nr and g^ ~ gN R indicate that the number of 
units (correlation blobs) required to form a topological 
constraint is on the order of ~ N^. 

The optimum free energy for large N(^> g\) is 
F I (ksT) ~ N/g^ , i.e., an order of the magnitude of ther- 
mal energy per topological strand. This translates into 
the osmotic pressure n top ~ FNr/R 3 ~ ksTcj)/ {a^g^) ~ 
(N R )~ 2 k B T4> 9 / 4 /a 3 . Adding this to the usual excluded- 
volume contribution Himear — ksTcj) 9 ^ 4 /a 3 [ill, the total 
osmotic pressure thus can be written as 

n = (i + N R 2 )u lmear (ii) 

where one may identify the small factor as the increase 
of the local excluded- volume effect, i.e., N^ 2 Sv/v in 
accordance to the des Cloizeaux conjecture [23j . 

To summarize, the present attempt strongly suggests 
the applicability of the excluded-volume analogy to the 
problem of dense ring solutions. It sets up a way to han- 
dle immense topological constraints, thus allowing one to 
capture the essential static properties in such systems. In 



particular, it provides us with a fairly accurate descrip- 
tion of the ring dimension over the range of semidilute 
to concentrated solution regimes. Other key quantities 
such as <fi s , Nr as well as the estimated onset length N* 
of the compact statics are also in line with reported nu- 
merical observations [9j-[l4j. It should be noticed that 
our starting point is similar in spirit to that in C-D the- 
ory, i.e., the identification of two competing contribu- 
tions in eq. (|2|). Rather, it is the difference in the basic 
structure of the free energy that leads to deep insights 
into the hierarchical spatial structure of ring solutions, 
hence a bridge to the CG concept. Further careful stud- 
ies are awaited to examine the validity and the limitation 
of the present phcnomenological approach as well as the 
usefulness of the concept such as the topological volume 
fraction, etc. As stated in ref. [l3|, there are various sit- 
uations in molecular systems, i.e., the collapse of a gel, 
the existence of chromosome territories, the compatibility 
enhancement Q etc. where topological constraints mat- 
ter. We hope that the present analysis provides valuable 
insights into such cases, too. T.S. thanks H. Nakanishi 
and A. Yoshimori for useful comments and discussions. 
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